Regge residues from DGLAP evolution 
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We show that combining forward and backward evolution allows to extract the residues of the 
triple-pole pomeron and of the other singularities for 10 GeV 2 < Q 2 < 1000 GeV 2 . In this approach, 
the essential singularity generated by the DGLAP evolution is considered as a numerical approxi- 
mation to a triple-pole pomeron. Using an analytical expression for the form factors, we reproduce 
the experimental data with a x 2 /dof of 1.02. This proves the compatibility between Regge theory 
and DGLAP evolution. The method used here enables us to evaluate the uncertainties on the gluon 
distribution which prove to be large at small x and small Q 2 . 



PACS numbers: 11.55.-m, 13.60.-r 
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I. INTRODUCTION 

Deep Inelastic Scattering (DIS) is a very important 
tool to study the fundamental behaviour of strong inter- 
actions. The numerous available data, especially the last 
ones from HERA, which have very small error bars, allow 
for a direct test of theoretical models. 

The work presented here puts together two major 
theoretical approaches: perturbative QCD (pQCD), 
through the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution and analytical S'-matrix theory, 
through Regge theory, on the other. Models using one 
or the other are numerous. On the one hand, there are 
a lot of well-known DGLAP fits to experimental data 
00000 which reproduce very well the experimen- 
tal data. However, in the usual parton distribution sets, 
each initial parton distribution presents its own singular- 
ities. As an example, in the MRST2001 parametrization 
0, we have 

xq(x, Ql) = A{\ + B\fx + Cx)(l - x) n "x e " , 

with e SGa = —0.26, e g — —0.33. Such a behaviour is 
not consistent with Regge theory. On the other hand, 
Regge theory establishes strong constraints on the large 
s (small x) behaviour of total cross sections, and Regge- 
compatible fits to experimental data are also numerous 
000 HH- In those fits, the residues of the Regge 
poles are unknown functions of Q 2 , and one can only 
predict F% while quarks and gluon distribution functions 
remain out of reach. 

In a previous paper we have shown that one can 
use a triple-pole pomeron to parametrize the initial par- 
ton distribution functions and let them evolve according 
to the DGLAP equation. In this approach, both quarks 
and gluons have the same singularity structure. Our 
model is therefore compatible both with Regge theory 
and with pQCD. In addition, we predict a gluon distribu- 
tion of the same order of magnitude as the one obtained 
in the usual DGLAP fits. 



We now want to explore whether it is possible to avoid 
the essential singularity coming from DGLAP evolution 
and extend the Regge fit to large Q 2 values. In other 
words, can we use the DGLAP evolution equations to 
find the Q 2 dependence of the Regge form factors ? To 
answer that question, the main obstacle lies in the fact 
that the DGLAP evolution generates an essential singu- 
larity which seems to be incompatible with the Regge 
ones. To solve that problem, we need to remember that 
the DGLAP evolution is only an approximation and that 
only a full resummation of the perturbative scries - in- 
cluding all order corrections to DGLAP, small-x and 
higher-twist effects - should give the correct analytical 
behaviour. From the large number of DGLAP fits, one 
expects corrections to be small and DGLAP to be a good 
approximation. We will therefore choose a Regge com- 
patible singularity structure and an initial scale Qq for 
DGLAP evolution, and then adjust the initial distribu- 
tion in order to reproduce the experimental data as well 
as possible. In this way, we obtain the residues at scale 
Qg and, repeating this procedure for different values of 
<5q we get the Regge form factors for the full range of 
the initial scale. In this method, it is important to note 
that we keep a common singularity structure, namely a 
triple-pole pomeron and an /, a2-reggeon, in all parton 
distributions. 



'Electronic address: g.soyez@ulg.ac.be 



The structure of the paper is as follow: in section 2, 
we introduce the basics of DGLAP evolution and Regge 
theory, we emphasise the tools needed in the rest of the 
paper and we summarise previous work. In section 3, 
we give the form of the initial distribution that we will 
use as initial condition for DGLAP evolution and that is 
obtained from Regge theory. Section 4 tells how we use 
DGLAP evolution to obtain the form factors and section 
5 gives the results of the fit. We conclude in section 6. 
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II. PERTURB ATIVE QCD AND REGGE 
THEORY 

A. Perturbative QCD 

In pQCD, the high-Q 2 behaviour of DIS is given by the 
DGLAP evolution equations [T^ . These equations intro- 
duce the parton distribution functions Qi(x, Q 2 ), q%(x, Q 2 ) 
and g(x, Q 2 ), which can be interpreted as the probability 
of finding, in the proton, respectively a quark, an anti- 
quark or a gluon with virtuality less than Q 2 and with 
longitudinal momentum fraction x. When Q 2 — ► oo, the 
Q 2 evolution of these densities (at fixed x) are given by 
the DGLAP equations 
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At leading order, P xy does not depend on Q 2 . The con- 
nection to F 2 is given by 

F 2 (x, Q 2 ) = x^2 4 [©(a;, Q 2 ) + Qi(x, Q 2 )] , (2) 

i 

where the sum runs over all quark flavours. 

The usual way to use this equation is to choose a set 
of initial distributions qi(x, Qq, a) to compute q%(x, Q 2 , a) 
using and to adjust the parameters a in order to re- 
produce the experimental data. This approach have al- 
ready been successfully applied many times 0, 0, 0, 0, Q 
and is often considered a very good test of pQCD. Nev- 
ertheless, these studies do not care about the singularity 
structure of the initial distributions, ending up with re- 
sults that disagree with Regge theory, and presumably 
with QCD. 



B. Regge theory 

Beside the predictions of pQCD, we can study DIS 
through its analytical properties. In Regge theory 
0,0, we consider amplitudes A(j,t) in the complex 
angular momentum space by performing a Sommcrfcld- 
Watson transform. In that formalism, the large-s de- 
pendence of amplitudes is determined by their singular- 
ity structure in j-plane. The singularities are functions 
of t only and this technique can be applied to the do- 
main cos(#i) 3> 1, where 6 t is the deviation angle in 
the center-of-mass frame for the t-channel amplitudes 
and where amplitudes are analytically continued to any 
cos((9 t ) value. For example, we can fit the DIS data or 
the photon structure function at large p.q = v (small x), 
and the total cross sections at large s. 

The models based on Regge theory 0, B use 
a pomeron term, reproducing the rise of the structure 



function (resp. cross sections) at small x (resp. large s), 
and reggeon contributions coming from the exchange of 
meson trajectories (a and /). 

We shall consider here the following parametrization 
for the pomeron contribution to F2 or to o~tat' 

a{Q 2 )\n 2 [v/v {Q 2 )\+c{Q 2 ), (3) 

corresponding to a triple pole in the complex j-plane 0] 

a 2aln(^ ) aln 2 (^ ) + c 



(7 - i) 3 (i - 1) 5 



i-i 



This seems to be the preferred phenomenological choice 
at Q 2 = 10]. Note that the upper expression, given 
in terms of v and Q 2 , can be rewritten in terms of Q 2 
and x = Q 2 /(2i/). So, the x behaviour of the structure 
functions is constrained by Regge theory. 
However, this approach has two problems: 

1. it only predicts F2, while a full set of parton distri- 
bution functions is needed, 

2. it does not give any information about the Q 2 de- 
pendence of the residues. 



III. MIXING DGLAP EVOLUTION WITH 
REGGE THEORY 

From a previous study [ill ], we have learned that the 
small-x behaviour of the parton distributions can be de- 
scribed by a triple-pole pomeron model at the initial scale 
Qq and then evolved using the DGLAP equation. We can 
summarise these results as follows: 

• It is possible to have the same singularity structure 
for quarks and gluons, in agreement with Regge the- 
ory. 

• One can use a triple-pole pomeron to parametrize 
the initial parton distributions. In particular, we 
predict a gluon distribution compatible with usual 
ones. 

• The DGLAP evolution generates an essential sin- 
gularity at j = 1. We expect that this singularity is 
purely of perturbative origin and that a full resum- 
mation of QCD should lead to a Regge-compatible 
singularity structure. 

• The F 2 residues were fixed to reproduce results 
from the soft QCD fit 0] which goes down to 
Q 2 = 0. 

We show here that we can use this trick to obtain form 
factors which are compatible with DGLAP evolution. We 
present in this section the model used to extract these 
residues as a function of Q 2 . In summary, always consid- 
ering that DGLAP evolution does not lead to a relevant 
singularity structure but only to an approximation of the 



3 



parton distributions, we consider that the initial distri- 
bution fitted by starting the evolution at the scale Qq 
gives the residues at Qq. We just need to repeat this by 
varying Qq to get the residues over the whole Q 2 range. 
Note that since we want to predict the F% residues, we do 
not constrain F2 with a soft QCD fit as it was done in 
However, it is known that the triple-pole pomeron 
model can be used to reproduce soft processes. There- 
fore, one knows that the residues found here can be ex- 
tended down to Q 2 = 0. Since we only want to study 
the domain where both DGLAP and Regge theory are 
expected to apply, we won't go into such an extension. 
Regge theory is expected to be valid for cos(# t ) » 1. 

In the case of DIS, cos(0 t ) = ^— 



where 

of the proton. The Regge region is thus 



is the mass 



>K, 
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FIG. 1: Number of points involved in the forward evolution 



with K a fixed number. We clearly see that the domain 
does not extends up to x = 1. However, in the DGLAP 
evolution equations the small- 2; domain is coupled to 
large- a: values and therefore we need a parametrization 
for the parton distribution functions at large x. The ideal 
case is certainly the introduction of powers of (1 — x), 
or of any polynomial vanishing at x = 1, in our initial 
distributions. This should take into account daughter 
trajectories in Regge theory and ensure that parton dis- 
tributions go to as x — > 1. Unfortunately, this has two 
drawbacks: 

1. it introduces a lot of additional parameters which 
are not directly related to the small- ir Regge be- 
haviour, 

2. a precise description of the parton distributions at 
large x requires more than simply fitting F% and 
splitting it into flavor singlet component, coupled 
to the gluons, and one additional non-singlet distri- 
bution. We need to introduce valence quark distri- 
butions and sea quark distributions. Fitting many 
kinds of experiments like F 2 P , F£, F 2 d , F£ N and 
xF£ N allows to constraint all these distributions 
separately. Again, introducing more parton distri- 
butions requires more parameters. Moreover, most 
of the experimental points except those from i^, 
lies outside the region of interest (0J . 

In order to keep the parameters of our model as closely 
linked as possible to the Regge domain of the evolution, 
we will use an external distribution for the large-x parton 
distributions. As done in we will use GRV98 at 
large x. We have already argued that the results do not 
depend much on that choice and that their leading singu- 
larity structure is not affected by the parametrization at 
large x. Practically, we shall use the GRV parametriza- 
tion for x > xgrv an d our analytical form at small x. 



IV. INITIAL DISTRIBUTIONS 

In order to reproduce F^ , we do not need the full set 
{<Zi>9ii3} of parton distributions. One can quickly check 
that the minimal set of distributions required is 

T = x [{u+ + c+ + t+) - (d + + s+ + &+)] , (5) 
E = x [(u+ + c++ t+) + (d+ + s+ + &+)] , (6) 
G = xg, (7) 

where q + = q + q for q = u,d, s, c, t, b. The T distribution 
is a flavor non-singlet distribution, its Q 2 evolution is 
decoupled from any other distribution and is given by 
a DGLAP equation with xP qq (x) as splitting function. 
The sea quark distribution S is flavor singlet and evolves 
coupled to the gluon distribution using the full splitting 
matrix. 

Using these distributions, the proton structure func- 
tion becomes 



F 2 = 



5E + 3T 
18 ' 



In order to perform a DGLAP evolution, we must find 
an analytical expression for the x dependence of these dis- 
tributions at a scale Qq, where evolution will start. Since 
we want to describe F 2 as a triple-pole pomeron and a 
reggeon (/, a 2 -trajectories), and since adding new singu- 
larities in the parton distributions disagrees with Regge 
theory, the T, S and G distributions will be described 
by: 

T(x,Ql) = a T log 2 (l/x) +6 T log(l/a;) + c T + d T x\ 
E(x,Ql) = a s log 2 (l/a;) + 6 s log(l/x)+c s +d s x' 7 , 
G(x,Ql) = a G \og 2 (l/x) + b G \og(l/x)+c G + d G x ri . 

Out of the 12 parameters in the above expressions, 
some may be removed. First of all, we expect that the 
pomeron, having vacuum quantum numbers, does not 
distinguish between quark flavours. This means that it 



4 



is decoupled from the T distribution and that we may set 
ax = bj- = ct = 0. Then, since the reggeon is mainly 
constituted of quarks, we shall neglect its contribution to 
the gluon distribution and set d G — 0. 

In addition, we ask our distributions to be continuous 
with GRV's at x = xqrv- This is used to fix aV, cs and 
cg- 

Finally, we have seen that we obtain much better re- 
sults if we multiply each distribution by some power of 
(1 — x). This new factor, effectively including daughter 
trajectories, is expected from Regge theory. As already 
pointed out, this factor is not directly connected to the 
small- x domain. One can see that its main effect is to 
reproduce an inflexion point present around x w 0.01 in 
the parton distributions. In order to minimise the num- 
ber of parameters, we used only two different powers in 
all distributions. We thus end with the following func- 
tions 

T(x,Ql) = d* T x 7 \l-x) b \ 

£(x,Q2) = as log 2 (l/ a; ) + 6 E log(l/x)+ C *(l-x)'' 1 

+ d^il-xf 2 
G(x,Ql) = a G \og 2 {l/x) + b G \og(\/x)+c G {l-x) h \ 

The parameters with a superscript * are constrained by 
continuity with GRV parametrisation and 7 free param- 
eters remain: as, &e ; ^s, ad b G , &i an d 62. 

These parameters must depend on Qq and we will see 
in the next section how we can find their Q 2 dependence. 



V. FITTING THE FORM FACTORS 

It is well-known that DGLAP evolution generates an 
essential singularity at j = 1. From that point of view, 
finding the Q 2 dependence of the residues of a triple- 
pole pomeron may seem impossible. However, DGLAP 
evolution always involves a finite order calculation of the 
splitting functions and, moreover, it resums large Q 2 cor- 
rections while resumming large s (small x) contributions 
should be preferable for a comparison with Regge the- 
ory. Consequently, we argue that the singularity struc- 
ture generated by DGLAP evolution is physically irrel- 
evant and that only a full resummation of QCD should 
end up on a correct Regge behaviour. 

One must note that if we consider a resumation a la 
BFKL, we have a stable square- root branchpoint instead 
of an essential singularity. Furthermore, even in this case 
we still need to unitarize the pomeron, which may lead 
to another singularity structure. 

So, which singularity do we need to take ? How can we 
use the DGLAP evolution to calculate residues in Regge 
theory ? Our point of view is that, although it doesn't 
provide a correct Regge behaviour, DGLAP evolution 
can still be considered as a good numerical approxima- 
tion. In such a case, the only place where Regge theory 
can be used is for the initial distribution. We end up 




FIG. 2: Form factors as functions of Qq 



with a Regge-compatible distribution at the initial scale 
Qq and a DGLAP approximation everywhere else. 

As we shall see, obtaining the Regge residues from this 
technique is a quite straightforward. 
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A. Scheme 1 

As initial condition for DGLAP evolution, we shall use 
the distributions obtained in section IIVI for x < xqrv 
and the GRV98 distribution for larger x values. The 
parameters are considered as functions of Q 2 and they 
can be found as follows: 

1. choose an initial scale Qq, 

2. choose a value for the parameters in the initial dis- 
tribution, 

3. compute the parton distributions at larger Q 2 using 
DGLAP evolution equations, 

4. repeat |21 and El until you find the best values of the 
parameters reproducing the experimental F 2 data 
for Q 2 > Qq and x < xqrv- 



5. 



This gives the residues at the scale Qq and we re- 
peat steps H to 0] in order to obtain the residues at 
all Q 2 values. 



This technique works quite well but has one drawback: 
calculation of the residues at the scale Qq relies on the 
experimental points satisfying 



Ql < Q 2 < Q\ 



,2 

max ■ 



Unfortunately, this domain depends on the initial scale 
and as we can see on Fig. ^ the number of points used 
for the fit decrease quite fast when Qq rises. 



B. Scheme 2 

There is a way to keep the same set of points in each 
fits: adding backward evolution. The previous scheme 
thus becomes 

1. choose an initial scale Qq, 

2. choose a value for the parameters in the initial dis- 
tribution, 



3. 



4. 



■5. 



compute the parton distributions for Qq < Q < 
Q max using forward DGLAP evolution and for 
Qmin < Q 2 < Qq usin g backward DGLAP evo- 
lution, 



repeat |21 and until you find the value of the pa- 
rameters reproducing the F 2 data for Q 2 > Q 2 lin 
and x < £grv- 

This gives the residues at the scale Qq and we re- 
peat steps n to 0] in order to obtain the residues at 
all Q 2 values. 



This way, we take into account the experimental points 
with <3 min < Q 2 < <5 max whatever the value of Ql is. 
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FIG. 3: xV n of th e fit as a function of the initial scale for 
DGLAP evolution 



Experiment 


nb. pts. 


x 2 


X 2 /n 


BCDMS[16J 


5 


20.9 


4.185 


E665[1JJ 


10 


8.0 


0.797 


Hl[18,i9, 20] 


240 


219.2 


0.913 


NMC[21j 


11 


17.1 


1.553 


ZEUS [22, 23] 


294 


306.3 


1.042 


Total 


560 


571.5 


1.020 



TABLE I: Experimental points in the domain ® and results 
of the analytical triple-pole fit 



VI. RESULTS 

We have chosen to use the second scheme with both 
forward and backward evolution. We have fitted the data 
in the region 



10 < Q 2 < 1000 GeV 2 
x < zgrv = 0.15, 

9 t ) = y^! > 



(8) 



cos 



49 GcV^ 
2m'i 



This method ensures that Q 2 is large enough to apply 
DGLAP evolution, that cos(# t ) is large enough to use 
Regge theory and that the large- a; domain is excluded 
from the fit. The limit on cos(6> t ) has been taken from 
the pure Regge fit of 6j. In this region, the experimental 
points are as shown in Tabled 



A. Form factors 

When we apply the method explained in section IV Bl 
using the parametric distributions obtained in section llVl 
we obtain the form factors presented in Fig. and cal- 
culated from 10 to 1000 GeV 2 . The \ 2 P er data point 
obtained for each fit, shown in Fig. |21 is of order 1 and 
the form factors appear to be smooth functions of Q 2 . 

Finally, when choosing xgrv, wc do not want to take 
the high- a; experimental points for reasons explained pre- 
viously, but we want to set xgrv as high as possible 
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to maximize the number of experimental point as pos- 
sible, particularly in the large Q 2 domain where DGLAP 
evolution is expected to work better. We found that 
2-grv — 0-15 w as a good compromise between those con- 
straints. 



B. Parton Distributions 



Once we have fitted the form factors, it is interesting 
to see if the distributions we obtain at a scale Q 2 , are 
the same if we start evolution at Q 2 or if we start at 
Q\ and evolve until Q\. As we see in Fig. 01 for the 
sea-quark distribution, the difference is very small. This 
shows clearly that the DGLAP essential singularity can 
be considered as a numerical approximation to a triple- 
pole singularity, at least for quarks. 

For the gluon case, presented in Fig. things are dif- 
ferent. One can see that backward evolution tends to pro- 
duce lower gluon distribution at small x. These can even 
become negative if we evolve to lower Q 2 values. How- 
ever, the structure function only depends on the quark 
distributions which is stable as shown in Fig. 0] Thus, 
from the phcnomcnological point of view, all these distri- 
butions provide a correct description of the data and are 
perfectly acceptable. The differences in the gluon distri- 
bution obtained for different values of Qq must therefore 
be considered as uncertainties on the gluon distribution. 
One clearly sees from Fig. |31that even at Q 2 — 10 GeV 2 , 
the uncertainties are very large. These large errors on 
the gluon distribution at small x and small Q 2 can be of 
prime importance in the study of LHC physics. 
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Q 2 



4.78 
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Using these analytical expressions, one can compare 
their predictions with the experimental data. We obtain 
a x 2 of 1.02 per data point, in a very good agreement with 
experimental measurements. The results are presented in 
Table □ and in Figs. E]and|7| 

In the above equations, we only wanted to find some 
analytical expression reproducing the form factors ex- 
tracted numerically from the DGLAP fits. We have not 
taken care of their analytical properties in Q 2 . To obtain 
better expressions, one should find an analytical equa- 
tion for the form-factors Q 2 -evolution. Since we have 
seen that the evolved distributions, containing an essen- 
tial singularity, are numerically close to the triple-pole 
distributions, one may hope that subleading corrections 
to the DGLAP equations will stabilize a triple-pole dis- 
tribution. 



D. F 2 C and F L 



C. Analytical form factors 



As a final step, we have tried to find analytical expres- 
sions for the form factors produced by the fit. Actually, 
it is much easier to work with analytical expressions re- 
producing results of Fig. |21 than to refit everything for 
each Q 2 values. In the Q 2 -range studied here, we have 
found that the following expressions, plotted in Fig. |3 
reproduce very well the form factors extracted from the 
DGLAP fit. 



Although only F% is included in the fit, it is interesting 
to see if our predictions agree with the measurements of 
F£ HH2II23 and F L Actually, these quantities are 
directly sensitive to the gluon distribution and are there- 
fore a good test of our predictions. The charm structure 
function is given by |26j 

/ «9«.0 2 + 4m;;)/(i/{,0 2 ), 

J ax 
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x x 



FIG. 4: E distribution at various Q 2 . Different curves correspond to different values of Qq. 




FIG. 5: gluon distribution at various Q 2 . Different curves correspond to different values of Qq. 




FIG. 6: F2 fit using analytical approximation of the residues (low Q 2 values) . 
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with 



/(x,Q 2 



(4- li)x 2 (1 - x) 



'X 

.2 

2m 2 r 



i (l-x) + fj,x 2 (l-3x)- fi 2 



L = log 



2x/i 
" 1-x 
1 + v' 



1-v 
l + 2fi. 



We have adopted a value of 1.3 GeV for the charm quark 
mass and the predictions for F% obtained from our model 
are presented in Fig. [SJ We see that the data are quite 
well reproduced. 

For the case of the longitudinal structure function, we 
are sensitive both to quarks and gluons: 



F L = 



G„ 



q—u,d,s,c 



4a s (Q 2 ) f 1 d£ fx 



3?r 



where 

G a 



2e 2 



a s (Q 2 + 4m 2 ) 



g(^Q 2 +±m 2 c ) 



X\ (XX 

1 - - ) v — L 



with m u — m,d = m s = and m c = 1.3 GeV. The results 
obtained from our parton distributions set are shown in 
Fig. and show good agreement with the data. 



VII. CONCLUSION 

In summary, we have seen in this paper that DGLAP 
evolution can be used to extract the residues of a triple- 
pole pomeron model. We chose triple-pole distributions 
for both quarks and gluons at an initial scale Q\ and 
we evolved them using DGLAP evolution. A compari- 
son with the experimental data gives the residues of the 
triple-pole pomeron at the scale Q\. In order to keep the 
same data points when Q\ varies, we evolved the initial 
distribution both forward and backward between Q^ aia 
and Q^ ax . We have repeated that fit for various values 
of Qq resulting in smooth form factors for the triple-pole 
residues. 





FIG. 9: Predictions for the longitudinal structure function using our parton distributions 
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In this model, quarks and gluons have the same sin- 
gularity structure, as suggested by Regge theory. The 
essential singularity generated by DGLAP evolution ap- 
pears to be a good numerical approximation to a triple 
pole for the F2 structure function. For the gluons, we 
found that there were quite large uncertainties in the 
gluon distribution at small x and small Q 2 . These un- 
certainties should certainly be taken into account in the 
LHC phenomenology. 

Finally, we were able to find analytical expressions for 
our form factors. A comparison of these expressions 
with experimental data gives a \ 2 P er point of 1.02, 
which proves the very good agreement between pQCD 
and Regge theory. We have also checked that our gluon 
distribution was acceptable by comparing our predictions 
for F% and Fl with the HI and ZEUS measurements. 

Fitting the parton distributions up to x = 1 in- 
volves many complications, therefore we used the GRV98 
parametrization for the large-x behaviour. As pointed 
out previously, this does not influence the small- x dis- 



tributions much. Extension of this model up to x — 1 
is left for future work. Since we know that Regge mod- 
els can reproduce soft-processes measurements, it should 
also be very interesting to check that we can also extend 
the form factors found in this paper down to Q 2 = 0. 

Finally, we used DGLAP evolution to extract the 
residues from the data. As we have seen, these are 
smooth functions of Q 2 and give numerical predictions 
close to the DGLAP essential singularity. It should there- 
fore be of prime interest to find an equation which sta- 
bilizes the triple-pole singularity and keeps the correct 
large-Q 2 behaviour. 
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